From Anderson to anomalous localization in cold atomic gases with effective 

spin-orbit coupling 
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We study the dynamics of a spin-orbit coupled Schrodinger particle with two internal degrees of 
freedom moving in a ID random potential. We show that this model can be implemented using 
cold atoms interacting with external lasers. Numerical calculations of the density of states reveals 
the emergence of a Dyson-like singularity at zero energy when the system approaches the quasi- 
relativistic limit of the random-mass Dirac model for large spin-orbit coupling. Simulations of the 
expansion of an initally localized wave-packet show a crossover from an exponential (Anderson) 
localization to an anomalous power-law behavior reminiscent of the zero-energy (mid-gap) state of 
the random-mass Dirac model. We discuss conditions under which the crossover is observable in an 
experiment and derive the zero-energy state, thus proving its existence under proper conditions. 
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Tightly confined ultracold atomic gases [l[ provide 
an ideal arena for studying a broad spectrum of one- 
dimensional (ID) quantum phenomena. Examples range 
from Hubbard-type models 0, ! pairing phenomena in 
Fermi gases Q to spin-orbit (SO) coupling for neutral 
atoms jp, El and quasi-relativistic physics Q, Q , gener- 
ated by the additional use of external laser fields. The 
same experimental techniques can also be used to address 
the effects of disorder in condensed matter systems. A 
prominent example is the demonstration of Anderson lo- 
calization of a Bose-Einstein Condensate (BEC) in a 
ID wave-guide (lp|.[l~lT|. One dimensional Anderson local- 
ization is caused by destructive interference in a weak, 
disordered potential (referred to as diagonal disorder) 
and leads to exponential localization of the particles 
15j. This behaviour may change in half- filled disordered 
metals [r| , random spin-Peierls and spin- ladder systems 
[l7^ . or as recently predicted in photonic systems with 
electromagnetically induced transparency [18J. Delocal- 
ized zero-energy (mid-gap) states can emerge in these 
systems showing a power-law behavior for the correla- 
tions due to a Dyson singularity in the density of states 
[l9| . Such anomalous localization originates in the chi- 
ral symmetry of the corresponding ID Hamiltonian and 
can be realized in the system with off-diagonal disorder 
known as a random-mass Dirac model or the fluctuating 
gap model (FGM) (23423- 

In cold atom systems disorder is typically induced by 
a random potential and is thus of diagonal type. Here we 
show that the combination of a random potential with SO 
coupling induced by the motion in space dependent laser 
fields can lead to effective off-diagonal disorder. Light- 
induced SO coupling has been shown to lead to an effec- 
tive Dirac dynamics in By investigating the density 



FIG. 1: (Color online) Left: Experimental sketch of the 
system. The condensate is placed in a tight cigar-shaped trap 
and driven by three lasers with Rabi frequencies fii, Q2 and 
Q3. Right: Internal tripod-type linkage pattern of the atoms 
exhibiting two dark-states without any contribution of the 
excited state |0). 



of states of the corresponding disorder model we derive 
conditions under which power-law localization can be ob- 
served and argue that they are indeed connected to the 
emergence of a Dyson singularity in the density of states 
(DOS). We show by simulating the time evolution of an 
initially localized wave packet that increasing the SO cou- 
pling drives a crossover from exponential (Anderson) lo- 
calization to an anomalous power-law localization. We 
note that a model similar to ours has been considered in 
24j, however only for the case of diagonal disorder and 
neglecting the kinetic energy part. As will be shown the 
latter can be neglected only under special circumstances 
and drastically alters the behavior of the system leading 
to much richer physics. 

We consider an ensemble of atoms exposed to three 
laser fields in a tripod-type linkage pattern (2j| as de- 
picted in Fig. (fTJ). The atoms are characterized by a de- 
generate manifold of the three ground states |1), |2), |3), 
coupled to a common excited state |0) via correspond- 
ing control lasers of equal wave- number k. We assume 
two of the control lasers to be counter-propagating along 
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the x-axis with Rabi frequencies = Q sin 9e~' lKX / \/2 
and H2 = O,sm0e lKX /v2, whereas the third laser prop- 
agates along the y-axis with 5I3 = ft cosde~ lKy , and 



fl = y£Li |^i| 2 denotes the total Rabi frequency. 
Following [2^1 we can write the interaction Hamilto- 
nian (h = 1) of the resulting tripod scheme as Hq = 

— Ei=i ffij|0)(i|+/i.c.^. This Hamiltonian has two dark- 
states, = -^"^(e^l) - e^ ma; |2)) and \D 2 ) = 
jpfT lKy cos0(e irex |l) +e~ iKX \2)) -sin0|3), with zero en- 
ergy and decoupled from the excited state |0) . We express 
a general state of the system in the dark state manifold 
according to |x(r)) = Yh=i ^i( T )\ A(r)), where Vi(r) 
and ip2 (r) are the wave functions corresponding to the 
two degenerate dark states. An effective equation for the 
center-of-mass amplitudes ipi{ r ) 1S then given by Q 
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where p is the momentum operator, m the atomic mass, 
and 9(r) = (-01 (r), -02(r)) T a two-component vector. 
The gauge potential A, known as the Mead-Berry con- 
nection [26U27|. arises from the coordinate-dependence of 
the dark-states and is given by Ak, n = i(Dk(v)\'V\D n (r)), 
The external potential has matrix elements \k,n = 
(Dk(r)\V\D n (r)) , where the potential in the bare basis 
was assumed to be diagonal, i.e., V = Ya=i ^*( r )K) 
The scalar potential reads &k,n = E;=i Afc l jA; j „/2m. 

By applying an additional strong transverse trapping 
potential to freeze out the transverse degrees of freedom, 
Eq. ((T|) reduces to a ID equation with 2 x 2-matrices 
A = -ncos8cr x , V = Diag[Vi,Vicos 2 + V 3 sin 2 6>], 
and 3? = k 2 /2m Diag[sin 2 0, sm 2 (26*)/4], where cr,, i € 
{x,y, z}, are the Pauli matrices. Expanding the square 
in Eq. ([1]) and choosing the external potentials as 
Vi = A(x) - n 2 /2m, and V 3 = -A{x) - k 2 cos 2 /2m, 
where A(x) is a detuning we arrive at 

<9# / p 2 \ 

* ~dt = \ 2m~ + VdCOS0 Vx ' &x + A ( a; )^) ®. ( 2 ) 

with vo = n/m. Eq. ^ describes a massive parti- 
cle with SO coupling moving in a scalar potential A (a;). 
The corresponding dispersion relation for a constant po- 
tential A(x) = S can easily be calculated and is shown in 
Fig. ([2|) for the case of 6 — (red dashed line) and small 
S (blue solid line). By changing the relative inten- 
sity of the control lasers one can externally control the 
value of cos# enabling an easy access to tune the effec- 
tive SO coupling. In case of a smooth potential and large 
particle-momenta, i.e., (p 2 )/2m 3> vo(p x ) cos0, the SO 
coupling becomes negligible and the problem reduces to 
two uncoupled massive Schrodinger particles moving in 
an external potential. In the opposite limit, i.e., for 
(p 2 )/2m <C vd(Px) cos8, Eq. ^ reduces to an effective 




FIG. 2: (Color online). Dispersion relation in the Schrodinger 
limit, with the inset representing the Dirac limit. The red 
dotted lines correspond to the case of vanishing potential, i.e. 
5 = 0, whereas the blue solid lines correspond to S = mwfj/8. 
The mixing angle is given by 9 = 0. 



Dirac equation for a particle with an effective speed of 
light c* = vd cos 8 and a smooth space-dependent mass 

A(x)/d 0,11. 

By applying Zeeman and Stark shifts, using, e.g., 
speckle potentials fiol [28| or incommensurate optical lat- 
tices a spatially fluctuating, random detuning, and 
thus random potential, can be generated. In the follow- 
ing we assume that the random potential is described by 
local, Gaussian white noise, 



A{x)A{x')=Tv D &{x-x'), A(x) = 0, 



(3) 



where the overbar denotes disorder average and all 
higher order correlation functions factorize. For van- 
ishing SO coupling one expects exponential localiza- 
tion of the par ticles according to the Anderson scenario 
0, ED, 0, IH, El]. On the other hand, any small SO 



coupling dominates in the region of small kinetic energy. 
Neglecting the kinetic energy term proportional to (p 2 ) 
in Eq. ([2]), rather than the SO coupling, one obtains 
an effective Dirac equation for a particle with a spatially 
random mass. This model, also known as the fluctuating 
gap model (FGM) [22, 23 1, is characterized by a Dyson 
singularity in the density of states (DOS) [30(, which is 
a consequence of the chiral symmetry of Eq. ^ without 
kinetic energy leading to a power-law localization. 

What determines the crossover from Anderson-like to 
anomalous localization? In the case of a random po- 
tential it is not immediately obvious in which parame- 
ter regimes Eq. ([2|) leads to Schrodinger or Dirac-like 
dynamics. To investigate this question we introduce di- 
mensionless units £ = Tx/vd and r = Tt resulting in 
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where A(£) = A/T. In the following we analyze the 
DOS in the system described by Eq. f4]). For a free 
Schrodinger particle the DOS has a singularity at u = 0. 
As can be seen from Fig. [3] the presence of the SO cou- 
pling shifts the band-edge of the spectrum away from 
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FIG. 3: (Color online). Density of states per unit length for a 
SO coupled massive Schrodinger particle in a random poten- 
tial with a finite disorder correlation length [33j for various 
values of the band edge mc*. As the band edge is moving 
closer to u = the washing out of the Dyson singularity of 
the free Dirac case (blue line) becomes more prominent, until 
the almost purely disordered Schrodinger case (magenta line) 
is reached. 
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Dirac limit mv^ — > oo the band edge moves to infin 
ity. In the pure Schrodinger case weak disorder leads 
to a smoothing of the band-edge peak [3l[, making the 
DOS analytic, which is associated with the emergence 
of exponentially localized states (Anderson localization) 
[32l |. To investigate the DOS for our SO-coupled sys- 
tem in the presence of disorder we numerically simulate 
Eq. ([4]). As will be the case in any experiment we as- 
sume for this a discretized model and furthermore a fi- 
nite disorder-correlation length (33|. (A smaller correla- 
tion length is accompanied by slower numerical conver- 
gence of the DOS [22J.) The mobility edge in momen- 
tum space associated with the finite disorder correlation 
length is however chosen large enough not to affect the 
results. Fig. [3] shows the DOS for different values of the 
dimensionless quantity Ao/mc 2 , where Ao characterizes 
the rms value of the Gaussian disorder amplitude in the 
discretized model. It is related to the continuum quan- 
tity via T = Ao/(v£>n kinks ), with n kinks being the impu- 
rity density. For mc 2 — > 00, i.e. in the Dirac limit, one 
recognizes a Dyson- like singularity at uj = 0. The emer- 
gence of a Dyson singularity can in general be taken as 
an indicator for anomalous, i.e. non-exponential, local- 
ization properties [l9| . In the random-mass Dirac model 
or FGM the Dyson singularity has been shown to lead 
to power-law correlations [2(| • As mc 2 ~ mv 2 D cos 2 9 de- 
creases the band-edge approaches the singularity. As a 
consequence the singularity is smoothend out. It should 
be noted at this point that a true singularity is present 
only in the exact Dirac limit, but as can be seen from Fig. 
(J3j) a pronounced peak survives as long as the SO cou- 
pling, i.e. mc 2 , is large. Thus one expects in this case a 
localization scenario where power-law correlations domi- 
nate for very long time scales. To quantify this we define 



following [22j the width Awb of the Dyson singularity 
as the minimum of the DOS, yielding Aujd = aT, where 
a = 0.6257.... Hence, to obtain dynamics associated with 
the Dyson singularity the condition 

mv 2 D cos 2 9 > T, (5) 

should be satisfied in order to ensure that the band-edge 
is sufficiently far away from the singularity. 

The emergence of the Dyson singularity affects the 
dynamics of the system around uj = 0. Experimen- 
tally this is reflected in a drastically different behav- 
ior of the density profile after a long-time expansion of 
an initially localized wave-packet. From Eq.Q we ex- 
pect that for cos 2 9 ^s> YJ mv 2 D the effective Dirac dy- 
namics along with the creation of anomalous power-law 
correlations dominates the system whereas in the op- 
posite limit we expect Anderson-like localization. The 
power-law behavior of correlations will also be reflected 
in the density distribution of an expanding wavepacket, 
which is experimentally much easier accessible. To see 
this we note, that in the power-law case there is no in- 
trinsic length scale. Thus we may assume that an ini- 
tially well localized wave-packet will show the same be- 
havior as a wavepacket with an initial delta-distribution 
>P(x, t = 0) = 8(x)x with x being a constant vector with 
unity length. The time-evolution of this wave-packet 



is given by *(x,i) = £ [<£*(0) • x] <t>n{x)e 



where 



4> n {x) is the stationary state of a particular disorder re- 
alization belonging to the state of energy u> n . It imme- 
diately follows that the long-time evolution is given by 
|*(x,i->oc)| 2 = ^^(O)! 2 !^^)! 2 . The right hand 

n 

side of this expression corresponds to the localization 
criterion defined in (l9| . and it is thus sufficient to in- 
vestigate the long-time behavior of the density in order 
to determine localization properties. In the Schrodinger 
limit all correlations of the type C n (x) = \<j> n (0)\ 2 \<j> n (x)\ 2 
decay exponentially [l^, i.e. C n (x) ~ e -\ x \/ L ^^ wnere 
the localization length L\ oc is at most of the order of the 
system size. In the FGM the typical localization length 
of correlation functions scales as L\ oc ~ |lne 1 2 [U, where 
e is the distance from the band-center. For small ener- 
gies this length is much larger than the system size and 
consequently the correlation functions should be simi- 
lar to the zero-energy correlation function decaying with 
a power law scaling as C n (x) ~ (r|a;|/ij£))~ 3 ^ 2 at large 
distances 2l|. Accordingly one expects quite different 



scaling behavior in the two different limits of large and 
small SO coupling. To confirm the predictions we numer- 
ically simulated the time-evolution according to Eq. ^ 
of an initially localized wave-packet of the form \P(x, t = 

0) = (2V2^L erf(L/V2L )) _1/2 cxp(-a; 2 /4L 2 ) (l,*) r 
of width Lq in a system of total length 2L. The results, 
shown in Fig. ([4}, depict the density profile for different 
SO coupling strengths and Y jmv 2 D =0.1 after a suffi- 
ciently long time evolution T . For cos 9 = (0.01, 0.05) one 
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FIG. 4: (Color online). Stationary density profiles of an ini- 
tially localized wave packet of size LoV/vd = 7.5 after a time 
evolution Ft = 80 according to eq. @ with T/mv% = 0.1 
after averaging over 100 independent realizations. The sys- 
tem size is 2LT/vn = 240 and flunks = 6-25F/vd, i.e. 
Ag/r 2 = 6.25. (Top left, right) Schrodinger limit with 
cos# = (0.01,0.05). The green (straight) line corresponds 
to an exponential fit showing the expected Anderson local- 
ization. The inset shows the same density in log-log. (Bot- 
tom left) Density in the crossover regime with cos# = 0.5 in 
log-log representation (semi-log in the inset) and a power-law 
fit (green straight line) with exponent 3/2. (Bottom right) 
Dirac limit with cos 9 = 1 in log-log; the power-law behavior 
is clearly visible for large distances £ 3> 1. 



can clearly see the exponentially localized wings of the 
density reminiscent of Anderson localization in cold atom 
systems [hJ EH, 0| . Increasing the SO coupling strength 
leads to an intermediate regime at cos 6* = 0.5 where nei- 
ther exponential nor power-law localization can be deter- 
mined. Finally for strong SO coupling, i.e. cos# = 1 a 
fit with a power-law reveals that the density correlation 
scales as the expected 3/2 power- law at large distances 
in agreement with the predictions of the FGM. 

To understand the 3/2 power- law exponent and show 
its connection to the formation of a zero-energy state, 
we perform a gauge transformation of Eq. with 
\P = cxp (— imc^xa x ) $ and subsequently substitute the 
Ansatz <J> = x+ cxp (imc*x) + X- exp (— imc„x) giving 
(details see [3J 
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dx A(x) 
dx 2 



+ f y ®a y )x, (6) 



where terms proportional to d^X 01 fast oscillating expo- 
nentials were dropped under the assumption T/mVp <C 
cos 2 8. Further we defined the four-component object 
X = (X+.X-) T and a,, n, i t {x,y,z}, act on the 
momentum or internal degree of freedom, respectively. 
Equation © is a generalization of the model considered 
in (2(i 21 1 . It possesses a zero-energy (mid-gap) state 



given by x( x ) ~ ex P \ ±a/2c* J dyA(y) > Xo- The ma- 

l -°° J 

trix a = (f y 55 a z + f x ® a y ) has eigenvalues (1, —1, 0, 0). 

Choosing xo to be an eigenvector with eigenvalue ±1 
we obtain a power-law intensity correlation scaling as 



2l| in full agreement with the 



c(x) ~ (r|x|/^)~ 3 / 2 [2C 

numerical simulations. 

In summary, we have investigated the one-dimensional 
dynamics of a SO-coupled massive Schrodinger particle 
subject to a (5-correlated disorder potential. For weak 
SO coupling the system is equivalent to two independent 
Schrodinger particles with diagonal disorder. In the op- 
posite limit the system is described by the random-mass 
Dirac model with off-diagonal disorder. The model can 
be implemented with current state-of-the-art techniques 
by using atomic dark-states in cold atom systems. We 
showed by calculating the systems' DOS and direct nu- 
merical simulation of the time evolution of an expanding 
wave packet that there is a crossover from an exponential 
(Anderson) localized regime to a power-law regime gov- 
erned by a Dyson singularity when varying the strength 
of the SO coupling relative to the disorder strength. The 
crossover can be observed by expansion of an initially lo- 
calized wavepacket with an appropriately chosen width 
of the momentum distribution in a random potential. 

Future investigations should include effects of interac- 
tions, which will allow for studying models such as the 
relativistic Thirring model [35| for fermions as well as for 
bosons with a random mass. Interactions may well accel- 
erate the crossover to non-exponential behavior as they 
tend to delocalize particles. 
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